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We investigate the distribution of energy density in a stationary self-reproducing 
inflationary universe. We show that the main fraction of volume of the universe in a 



;_, ■ state with a given density p at any given moment of time t in synchronous coordinates 

is concentrated near the centers of deep exponentially wide spherically symmetric holes 
in the density distribution. A possible interpretation of this result is that a typical 
observer should see himself living in the center of the world. Validity of this inter- 
pretation depends on the choice of measure in quantum cosmology. Our investigation 
suggests that unexpected (from the point of view of inflation) observational data, such 
as possible local deviations from fi = 1, or possible dependence of the Hubble constant 
on the length scale, may tell us something important about quantum cosmology and 
particle physics at nearly Planckian densities. 
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1. In this paper we are going to show that according to a very wide class of inflationary theories, 
the main fraction of volume of the universe in a state with a given density p at any given moment 
of time t (during or after inflation) should be concentrated near the centers of deep exponentially 
wide spherically symmetric holes in the density distribution. 

Observational implications of this results depend on its interpretation. If we assume that we 
live in a part of the universe which is typical, and by "typical" we mean those parts of the universe 
which have the greatest volume with other parameters (time and density) being equal, then our 
result implies that we should live near the center of one of the holes in the density distribution. 
There should be many such holes in the universe, but each of them should be exponentially wide. 
Therefore an observer living near the center of any such hole will see himself "in the center of the 
world" . 

One should clearly distinguish between the validity of our result and the validity of its inter- 
pretation suggested above. Even though the effect by itself is rather surprising (its existence was 
first conjectured in |[^), we think that we can confirm it by several different methods. Meanwhile 
the validity of its interpretation is much less clear. Until the interpretation problem is resolved 
in one way or another, there will remain an open possibility that inflationary cosmology predicts 
that we should live in a center of a spherically symmetric well. As we will see, this possibility is 
closely related to the question whether or not fi = 1 in inflationary universe. 

2. It is well known that inflation makes the universe locally flat and homogeneous. On the other 
hand, during each time interval if~^ = ./ „ y,,-. the classical scalar field experiences quantum 
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jumps of a typical amplitude 6(j) ~ 2^- (Here H is the Hubble constant; we use the system of 
units Mp = 1.) These jumps lead to small density perturbations with almost flat spectrum 0. 
If we follow any particular inflationary domain, we can predict its behavior and calculate the 
typical amplitude of density perturbations in it. 

However, this program can be accomplished by the methods of ref. only if the value of 
the scalar field in the initial domain is small. Meanwhile if one starts with a domain containing 
a sufficiently large field (p, quantum fiuctuations 6(f) become so large that they may give rise to 
an eternal process of self-reproduction of new infiationary domains. For example, in the theory 
^0^ the self-reproduction occurs in domains containing the field > 0* ~ \~^^^ 0. Such 
a domain very soon becomes divided into exponentially many domains containing all possible 
values of the scalar field (i.e. all possible values of its energy density) at the same moment of 
tim,e. The knowledge of initial conditions in the original domain does not allow one to predict 
density perturbations or even mean density inside each of these new domains. Instead of that, 
one may try to study those domains which may look typical (i.e. most abundant) in the context 
of the new cosmological paradigm. 

In particular, one may ask the following question. Suppose that we have one infiationary 
domain of initial size H~^, containing scalar field (p > (p*. Let us wait 15 billion years (in 
synchronous time t in each part of this domain) and see, what are the typical properties of those 



parts of our original domain which at the present moment have density 10 ^^g-cm ^. As we will 
show, the answer to this question proves to be rather unexpected. 

This domain exponentially expands, and becomes divided into many new domains of size 
H~^, which evolve independently of each other. In many new domains the scalar field decreases 
because of classical rolling and quantum fluctuations. The rate of expansion of these domains 
rapidly decreases, and they give a relatively small contribution to the total volume of those parts 
of the universe which will have density 10~^^g-cm^'^ 15 billion years later. Meanwhile those 
domains where quantum jumps occur in the direction of growth of the field gradually push this 
field towards the upper bound where inflation can possibly exist, which is presumably close to 
the Planck boundary ^(0p) ~ 1. Such domains for a long time stay near the Planck boundary, 
and exponentially grow with the Planckian speed. Thus, the longer they stay near the Planck 
boundary, the greater contribution to the volume of the universe they give. 

Those domains which 15 billion years later evolve into the regions with density 10~^^g-cm~^ 
cannot stay near the Planck boundary for indefinitely long time. However, they will do their 
best if they stay there as long as it is possible. In fact they will do even better if they stay near 
the Planck boundary even longer, and then rush down with the speed exceeding the speed of 
classical rolling. This may happen if quantum fluctuations coherently add up to large quantum 
jumps towards small (j). This process in some sense is dual to the process of perpetual climbing 
up, which leads to the self-reproduction of inflationary universe. 

Of course, the probability of large quantum jumps down is exponentially suppressed. However, 
by staying longer near the Planck boundary inflationary domains get an exponentially large 
contribution to their volume. These two exponential factors compete with each other to give us 
an optimal trajectory by which the scalar field rushes down in those domains which eventually 
give the leading contribution to the volume of the regions of a given density p at a given time 
t. From what we are saying it should be clear that the quantum jumps of the scalar field along 
such optimal trajectories should have a greater amplitude than its naively estimated value ^, 
and they should preferably occur in the downwards direction. As a result, the energy density 
along these optimal trajectories will be smaller than the energy density of their lazy neighbors 
which prefer to slide down without too much of jumping. This creates holes in the distribution of 
energy density. We are going to show that at any given time t most of the volume of the regions 
of the universe containing matter of any given density p lies very close to the centers of these 
holes. 

3. The best way to examine this scenario is to investigate the probability distribution Pp(0, t) 
to find a domain of a given physical volume (that is why we use notation Pp) in a state with a 
given field at some moment of time t. The probability distribution Pp{(f), t) obeys the following 
diffusion equation (see [|T| and the references therein): 

dt 2 90 V 27r (90 V 27r ^^ 3/7(0) pJ^ ^'^^ p ^ > 



Note, that this equation is vahd only during inflation, which typically occurs within some limited 
interval of values of the field 0: (pmin < 4> < 4>max- In the simplest versions of chaotic inflation 
model 4>min ~ 1, and 4>max is perhaps close to the Planck boundary 0p, where V{(f)p) = 1. To find 
solutions of this equation one must specify boundary conditions. Behavior of solutions typically 
is not very sensitive to the boundary conditions at the end of inflation at (pmin, it is sufficient to 
assume that the diffusion coefficient (and, correspondingly, the double derivative term in the r.h.s. 
of equation (|l])) vanishes for < (pmin- The conditions near the Planck boundary = 0p play a 
more important role. In this paper we will assume, following [||], that there can be no inflation at 
V{(f>) > 1, and therefore we will impose on Pp the simplest boundary condition -Pp(0, t)|<^>0p = 0. 
At the end of the paper we will discuss possible modifications of our results if other boundary 
conditions should be imposed at the upper boundary of the region of values of the field where 
infiation is possible. 

One may try to obtain solutions of equation (P in the form of the series Pp{4>, t) = J^'^i e^°* T!'s{(p)- 
In the limit of large time t only the term with the largest eigenvalue Ai survives, Pp{(j),t) = 
e^^*7ri(0). The function tti in the limit t —>■ oo will have a meaning of a normalized time- 
independent probability distribution to find a given field in a unit physical volume, whereas the 
function e''*^* shows the overall growth of the volume of all parts of the universe, which does not 
depend on in the limit t ^ oo. In this limit one can write eq. (|l|) in the form 

1 ^ (i0Ma_ (H!!M^^^,M , iL f ZM ^,A ^ ,n(,) . ..„) ^ a. ..„) . (2, 
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In this paper we will concentrate on the simplest theory with V{(j)) = |0^, H = w ^0^. Eq. 
for this theory reads |1|: 
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We have solved this equation both analytically and numerically, and have found that the eigen- 
value Ai is given by d{X) Hmax, where d{X) is the fractal dimension which approaches 3 in the limit 
A — > 0, and H^ax is the maximal value of the Hubble constant in the allowed interval of 0. In our 
case Hmax = H{(f)p) = 2w ^. Thus, in the small A limit one has Ai = 3if (0p) = 2\/Qtt ^ 8.68. 

Note, that the distribution tti depends on very sharply. For example, one can easily check 
that at small (at < A~^/^Ai ) the leading terms in eq. (^ are the second and the last onesQ. 
Therefore the solution of equation (0) at < A^^/^A^ is given by 

TTi ~0V— -^1 . (4) 

This is an extremely strong dependence. For example, tti ~ 0^-^'^° for the realistic value 
A ~ 10^^'^. All surprising results we are going to obtain are rooted in this effect. One of the 



"'The appearance of the new scale (/> ~ A ^/^ in the theory X(j)'^ is explained in [Q, and was also pointed out by 
Mukhanov (private communication). 



consequences is the distribution of energy density p. For example, during inflation p ^ j4>^- This 
implies that the distribution of domains of density p is Pp{p) ~ p^'^° . Thus at each moment of 
time t the self-reproducing inflationary universe consists of indefinitely large number of domains 
containing matter with all possible values of density, the total volume of all domains with density 



2p being approximately 10 times greater than the total volume of all domains with density p\ 
If this result is understood (see its discussion in [^), all the rest should be easy... 

Let us consider now all parts of inflationary universe which contain a given field at a given 
moment of time t. One may wonder, what was the value of this field in those domains at the 
moment t — H~^ ? The answer is simple: One should add to (j) the value of its classical drift 
(J)H~^. One should also add the amplitude of a quantum jump A0. The typical jump is given by 
(50 = ±^. At the last stages of inflation this quantity is by many orders of magnitude smaller 
than (j)H~^. But in which sense jumps ±^ are typical? If we consider any particular initial 
value of the field 0, then the typical jump from this point is indeed given by ±^. However, if 
we are considering all domains with a given and trying to find all those domains from which 
the field could originate back in time, the answer may be quite different. Indeed, the total 
volume of all domains with a given field at any moment of time t strongly depends on 0: 



Pp{4>) ~ vri(0) ~ 0V ~ 1 ^ (f)^^ ^ see eq. (^). This means that the total volume of all domains 
which could jump towards the given field from the value + A0 will be enhanced by a large 
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additional factor p ij,\ ~ (-'-"'" ~^) ^ • ^^ ^^^ other hand, the probability of large jumps 
A0 is suppressed by the Gaussian factor expf — ^"^^^2*^ )• One can easily verify that the product 
of these two factors has a sharp maximum at A0 = Ai0 ■ ^, and the width of this maximum is 
of the order ^. In other words, most of the domains of a given field are formed due to jumps 
which are greater than the "typical" ones by a factor Ai0 ± 1. 

4. This result is very unusual. We became partially satisfied by our understanding of this effect 
only after we confirmed its existence by four different methods, including computer simulations. 
A detailed description of our results will be contained in a separate publication [^. Here we will 
briefly describe the computer simulations which we have performed. 

We have studied a set of domains of initial size H^^ filled with large homogeneous field 0. 
We considered large initial values of 0, where the self-reproduction of inflationary domains is 
possible. From the point of view of stochastic processes which we study, each domain can be 
modelled by a single point with the field in it. Our purpose was to study the typical amplitude 
of quantum jumps of the scalar field in those domains which reached some value 0o close to 
the end of inflation. 

Each step of our calculations corresponds to a time change At = uHq^. Here Hq = if (0o), 
and u is some number, u < 1. The results do not depend on u if it is small enough. The 
evolution of the field in each domain consists of several independent parts. First of all, the field 
evolves according to classical equations of motion during inflation. Secondly, it makes quantum 



jumps by 5(j) = ^J-^jf- sinrj. Here Tj is a set of random numbers, which are different for each 
inflationary domain. To account for the growth of physical volume of each domain we used 
the following procedure. We followed each domain until its radius grows two times, and after 
that we considered it as 8 independent domains. If we would continue doing so for a long time, 
the number of such domains (and our distribution Pp{(j),t)) would become exponentially large. 
We will describe in the procedure which we used to overcome this difficulty. In accordance 
with our condition Pp(0, t)|<^>0p = 0, we removed all domains where the field jumped to the 
super-Planckian densities ^(0) > 1. After a sufficiently large time t the distribution of domains 
followed by the computer with a good accuracy approached the stationary distribution tti which 
we have obtained in [Q by a completely different method. We used it as a consistency check for 
our calculations. 

We kept in the computer memory information about all jumps of all domains during the 
last time interval H^^. This made it possible to evaluate an average sum of all jumps of those 
domains in which the scalar field became smaller than some value within the last time interval 
Hq^. Naively, one could expect this value to be smaller than ■^, since the average amplitude 
of the jump is ■^, but they occur both in the positive and negative directions. However, our 
simulations confirmed our analytical result A0 = Ai0 ■ ■^. In other words, we have found that 
most of the domains which reach the hypersurface = 0o within a time interval At = Hq^ do it 
by rolling accompanied by persistent jumps down, which have a combined amplitude Ai0o times 
greater than ■^. 

Even though we obtained this result in the context of the theory ^0^, it can be easily gener- 
alized. One can show that for any infiationary theory in the regime where one can neglect the 
first (diffusion) term in the l.h.s. of equation (|^) (which is usually possible at the last stages of 
infiation), most of the domains of a given field are formed due to jumps which are greater than 
the "typical" ones by a factor ^lyrWr [|]- 

5. As we already mentioned, the probability of large fiuctuations should be suppressed by the 
factor expf — ^^2 ); which in our case gives the suppression factor ~ exp(— 10^). It is well known 
that exponentially suppressed perturbations typically give rise to spherically symmetric bubbles 
0. Note also, that the Gaussian distribution suppressing the amplitude of the perturbations 
refers to the amplitude of a perturbation in its maximum. Let us show now that the main part 
of the volume of the universe in a state with a given (or with a given density p) corresponds to 
the centers of these bubbles. 

Consider again the part of the universe with a given at a given time t. We have found that 
most of the jumps producing this field during the previous time interval H~^ occurred from a 
very narrow region + -|: + Ai0^ of values of the scalar field. The width of this region was found 
to be of the order of ^, which is much smaller than the typical depth of our bubble A0 ~ -^i0^- 
Now suppose that the domain containing the field appears not at the center of the bubble, but 
at its wall. This would mean that the field near the center of the bubble is somewhat smaller than 
0. Such a configuration could be created by a jump from (j) + -fj + -^i^^ o^ily if the amplitude 
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of the jump is greater than Ai0^. However, we have found that the main contribution to the 
volume of domains with a given is produced by jumps of an amphtude Ai0^ ± ^, the greater 
deviation from the typical amplitude Ai0^ being exponentially suppressed. This means that the 
scalar field can differ from its value at the center of the bubble by no more than the usual 
amplitude of scalar field perturbations ^, which is smaller than the depth of the bubble by a 
factor (Ai0)~^. Thus, the main fraction of the volume of the universe with a given (or with a 
given density of matter) can be only slightly outside the center, which may lead to a small dipole 
anisotropy of the microwave background radiation. 

We should emphasize that our results are based on the investigation of the global structure 
of the universe rather than of the structure of each particular bubble. If one neglects that the 
universe is a fractal and looks only at one particular bubble (i.e. at the one in which we live 
now), then one can easily see that inside each bubble there is a plenty of space far away from 
its center. Therefore one could conclude that there is nothing special about the centers of the 
bubbles. However, when determining the fraction of domains near the centers we were comparing 
the volumes of all regions of equal density at equal time. Meanwhile, the density p^aii of matter 
on the walls of a bubble is greater than the density Pcenter in its center. As we have emphasized 
in the discussion after eq. (|), the total volume of all domains of density p„aii is greater than the 
total volume of all domains of density Pcenter by the factor (p„aii/pcenter)'^'^° • Thus, the volume 
of space outside the center of the bubble is not smaller than the space near the center. However, 
going outside the center brings us to the region of a different density, pwaii > Pcenter- Our results 
imply that one can find much more space with p = p„aii not at the walls of our bubble, but near 
the centers of other bubbles. 

6. The nonperturbative jumps down should occur on all scales independently. At the earlier 
stage, when the time changed from t — 2H^^ to t — H^^, the leading contribution to the volume 
of the universe of a given density should be also given by domains which are jumping down by 
A0 ~ Ai0^. This implies that on each new length scale different from the previous one by 
the factor of e (one new e-folding) our trajectory runs down creating a local depression of the 
scalar field equal to Ai0 ^. One may visualize the resulting distribution of the scalar field in the 
following way. At some scale r the deviation of the field from homogeneity can be approximately 
represented as a hole of a radius r with the depth Ai0^. Near the bottom of this hole there 
is another hole of a smaller radius e~^r and approximately of the same depth Ai0^. Near the 
center of this hole there is another hole of a radius e~^r, etc. Of course, this is just a discrete 
model. The shape of the smooth distribution of the scalar field is determined by the equation 

A,0^ = ./AA,0^ (5) 



dlnrH ^ 2tx 
which gives 

02(r) ^ 0!(2) for r> H-K (6) 

l-Ai02(O)y^lnri7 

Note that 0(r) ^ 0(0) for r < H^^ (there are no perturbations of the classical field on this scale). 



This distribution is altered by the usual small perturbations of the scalar field. These per- 
turbations are responsible for the large nonperturbative jumps down when they add up near the 
center of the hole, but at a distance much greater than their wavelength from the center of the 
hole these perturbations have the usual magnitude ^. Thus, our results do not lead to consider- 
able modifications of the usual density perturbations which lead to galaxy formation. However, 
the presence of the deep hole (P) can significantly change the local geometry of the universe. 



In the infiationary scenario with V{(j)) = jcj)'^ fiuctuations which presently have the scale 
comparable with the horizon radius r^ ~ 10^^ cm have been formed at ~ 5 (in the units 
Mp = 1). As we have mentioned already Ai ~ 2-\/67r ~ 8.68 for our choice of boundary conditions 
[]I]]. This means that the jump down on the scale of the present horizon should be Ai0 ~ 40 
times greater than the standard jump. In the theory jcp^ the standard jumps lead to density 
perturbations of the amplitude — ~ 2v|A7r ^3 ^ g . xq-5 j^j^ ^-j^g normalization of ||7]]). Thus, 
according to our analysis, the nonperturbative decrease of density on each length scale different 
from the previous one by the factor e should be about — ~ Xi ^^'^ cp^ ~ 2 ■ 10^^. This allows one 
to evaluate the shape of the resulting hole in the density distribution as a function of the distance 
from its center. One can write the following equation for the scale dependence of density: 

1 dp 2a/6A^ 4 

= -Ai — , (7) 
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ro 



where r is the distance from the center of the hole. Note that = ^ (ln^)2 in the theory 

^0^ 0. Here tq corresponds to the smallest scale at which infiationary perturbations have been 
produced. This scale is model-dependent, but typically at present it is about 1 cm. This yields 



Ap ^ p{r) - p{ro) ^ 2AiV2A ^^3 r_ 
Pc ~ p{ro) SttVStt ro ' 

This gives the typical deviation of the density on the scale of the horizon (where In ^ ~ 60) from 
the density at the center: ^ ~ 750 ■ ^ ~ 4 ■ lO"^. 

7. It is very tempting to interpret this effect in such a way that the universe around us becomes 
locally open, with 1 — fi ~ 10~^. Indeed, our effect is very similar to the one discussed in 
1^, 0, where it was shown that the universe becomes open if it is contained in the interior of a 
bubble created by the 0(4) symmetric tunneling. Our nonperturbative jumps look very similar 
to tunneling with the bubble formation. However, unlike in the case considered in |^, 0, our 
bubbles appear on all length scales. This removes the problem of fine-tuning, which plagues the 
one-bubble models of the universe ^. However, it rises another problem, which may be even 
more complicated. 

Note, that the results discussed above refer to the density distribution at the moment when 
the corresponding wavelengths were entering horizon. At the later stages gravitational instability 
should lead to growth of the corresponding density perturbations. Indeed, we know that density 



perturbations on the galaxy scale have grown more than 10^ times in the linear growth regime 
until they reached the amplitude -^ ~ 1, and then continued growing even further. The same can 
be expected in our case, but even in a more dramatic way since our "density perturbations" on 
all scales are about 40 times greater than the usual density perturbations which are responsible 
for galaxy formation. This would make the center of the hole very deep; its density should 
be many orders of magnitude smaller than the density of the universe on the scale of horizon. 
Moreover, the center would be devoid of any structures necessary for the existence of our life. 
Indeed, on each particular scale the jump down completely overwhelms the amplitude of usual 
density perturbations. The bubble cannot contain any galaxies at the distance from the center 
comparable with the galaxy scale, it cannot contain any clusters at the distance comparable with 
the size of a cluster, etc. 

This problem can be easily resolved. Indeed, our effect (but not the amplitude of the usual 
density perturbations) is proportional to Ai, which is determined by the maximal value of the 
Hubble constant compatible with inflation. If, for example, the maximal energy scale in quantum 
gravity or in string theory is given not by 10^^ GeV, but by 10^*^ GeV, then the parameter Ai will 



decrease by a factor 10 . Other ways to change Ai are described in [10|. Thus it is easy to make 



our effect very small without disturbing the standard predictions of inflationary cosmology. 

It is possible though that we will not have any problems even with a large Ai if we interpret our 
results more carefully. An implicit hypothesis behind our interpretation is that we are typical, 
and therefore we live and make observations in those parts of the universe where most other 
people do. One may argue that the total number of observers which can live in domains with 
given properties (e.g. in domains with a given density) should be proportional to the total volume 
of these domains at a given time. However, our existence is determined not by the local density 
of the universe but by the possibility for life to evolve for about 5 billion years on a planet of our 
type in a vicinity of a star of the type of the Sun. If, for example, we have density 10~^^g-cm~^ 
inside the center of the hole, and density 10~^^g-cm~^ on the horizon, then the age of our part of 
the universe (or, to be more accurate, the time after the end of inflation) will be determined by 
the average density 10~^^g-cm~^, and it will be smaller than 5 billion years. Moreover, as we have 
mentioned above, any structures such as galaxies or clusters cannot be formed near the centers 
of the holes. 

Thus, the naive idea that the number of observers is proportional to volume does not work 
at the distances from the centers which are much smaller than the present size of the horizon. 
Even though at any given moment of time most of the volume of the universe at the density 
10~2^g-cm~^ is concentrated near these holes, the corresponding parts of the universe are too 
young and do not have any structures necessary for our existence. Volume alone does not mean 
much. We live on the surface of the Earth even though the volume of empty space around us is 
incomparably greater. 

On the other hand, the naive relation between the number of observers and the volume in a 
state of a given density may work nicely for the matter distribution on the scale of the present 
horizon. These considerations may make the distribution (H) homogeneous near the center not up 
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to the scale ro ~ 1 cm, but up to some much greater scale r ^ 10^^ cm. The resulting distribution 
may appear relatively smooth. It may resemble an open universe with a scale-dependent effective 
parameter Q{r). 

An additional ambiguity in our interpretation appears due to the dependence of the distribu- 
tion Pp on the choice of time parametrization. Indeed, there are many different ways to define 
"time" in general relativity. If, for example, one measures time not by clock but by rulers and 
determines time by the degree of a local expansion of the universe, then in this "time" the rate 
of expansion of the universe does not depend on its density. As a result, our effect is absent in 
this time parametrization |l[]. In this paper we used the standard time parametrization which 
is most closely related to our own nature (time measured by number of oscillations rather than 
by the distance to the nearby galaxies). But maybe we should use another time parametriza- 
tion, or even integrate over all possible time parametrizations? The last possibility, suggested 
to us by Robert Wagoner, is very interesting, since after the integration we would obtain a 
parametrization-independent measure in quantum cosmology. Even more radical possibility is 
the integration over all possible values of time. However, the corresponding integral diverges in 



the limit t — > cxd |]T0[. Right now we still do not know what is the right way to go. We do not 
even know if it is right that we are typical and that we should live in domains of the greatest 
volume, see the discussion of this problem in ||l], |10 . 



Until the interpretation problem is resolved, it will remain unclear whether our result is just 
a mathematical curiosity similar to the twin paradox, or it can be considered as a real prediction 
of properties of our part of the universe. On the other hand, at present we cannot exclude the 
last possibility, and this by itself is a very unexpected conclusion. Few years ago we would say 
that the possibility that we live in a local "center of the world" definitely contradicts inflationary 
cosmology. Now we can only say that it is an open question to be studied both theoretically 
and experimentally. As we have seen, the relation of our result to the properties of our part of 
the universe depends on interpretation of quantum cosmology, and the amplitude of our effect 
is very sensitive to the properties of the universe at nearly Planckian densities. This suggests 
that unexpected (from the point of view of inflation) observational data, such as possible local 
deviations from fi = 1, or possible dependence of the Hubble constant on the length scale, may tell 
us something important about quantum cosmology and particle physics near the Planck density. 
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